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ABSTRACT 


The utility of using a numerical weather prediction (NWP) forecast model as an 
input to a simple ocean model for planning during reactive situations is studied. An 
oceanographic experiment called the Maud Rise Nonlinear Equation of State Study 
(MaudNESS) was conducted from June to September of 2005 over the Maud Rise in the 
eastern Weddell Sea. Archived Antarctic Mesoscale Prediction System (AMPS) Polar 
MM5 forecast fields from MaudNESS were compared to observed conditions during 
MaudNESS. AMPS was found to have problems with cloud and moisture parameters, 
but represented the overall synoptic situation. AMPS forecast and observed forcing 
fields (as well as increased values for both) were input into a simple one dimensional 
ocean model at three locations in the Maud Rise area of differing stability. The ocean 
model was found to have good utility as a planning tool for short term reactive situations 


where a high degree of accuracy is not needed. 
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I. INTRODUCTION AND BACKGROUND 


The purpose of this research was to explore the feasibility of inputting numerical 
weather prediction (NWP) model forecast data into a simple ocean model, and then using 
the results of the ocean model in a reactive situation (best placement of a research vessel 
for investigating a desired oceanographic event, or best placement of fleet assets during 


an anti-submarine warfare (ASW) prosecution). 


The opportunity to conduct this research came after the conclusion of the Maud 
Rise Nonlinear Equation of State Study (MaudNESS). MaudNESS took place during 
June to September 2005 (Southern Hemisphere winter season) in the eastern Weddell Sea 
(see Figure 1) in the vicinity of the Maud Rise (a seamount in the eastern Weddell Sea). 
The seamount rises from the seafloor (at a depth of approximately 5000 meters) to 
approximately 1600 meters below the surface of the ocean, and is approximately 100 
kilometers in diameter. Muench et al. (2001) describe an isolated column of water 
partially trapped over the Maud Rise (resulting from impingement of combined Weddell 
Sea Warm Deep Water (WDW) and Circumpolar Deep Water (CDW)) which is 
characterized as a Taylor column. The Taylor column is colder, less saline, and slightly 
denser than surrounding waters, and tends to block horizontal flow over the Maud Rise. 
They describe a "halo" around the Maud Rise (ring shaped closed circulation about the 
Rise of relatively warm and saline combined WDW and CDW), which is the result of the 
regional southwestward flow being accelerated around the periphery of the Taylor 
column. Surrounding regional flow then continues southwestward playing a role in the 


formation of a warm pool southwest of the Maud Rise. 


The area around the Maud Rise has been known to be the site of formation for 
large polynyas. The Weddell Polynya formed there in the mid 1970's and persisted for 
several seasons, with cooling effects on WDW lasting for at least a decade. Another 
polynya was seen in satellite imagery late in the winter season of 1994. WDW directly 


feeds Antarctic Deep Water (AADW) which is extremely important to global 


thermohaline structure, therefore potentially significantly altering global deep water 


production. 


It is believed that thermobaricity (dependence of thermal expansion coefficient of 
seawater on pressure), cabbeling (dependence of thermal expansion coefficient of 
seawater on temperature), and other nonlinearities in the equation of state for seawater 
play an important role as preconditioners for deep convection to occur. In the case of the 
area around the Maud Rise, this deep convection is sometimes strong enough to 
counteract the stabilizing effect of buoyancy flux from ice melt (positive salinity flux). 
McPhee states in the MaudNESS cruise report (obtained from http://fish.cims.nyu.edu 
/project_maudness/field_campaign.html) that: 

The Maud Rise Nonlinear Equation of State Study (MaudNESS) was 


developed to investigate upper ocean mixing in the low-stability regime 
found in the Maud Rise region. 


Stanton describes on his website (http://www.oc.nps.navy.mil/~stanton/thermo 
/Maudness/IntroductionMain.html) "the objectives of the experiment is to answer two 


sets of questions: 


1.) What are the mechanisms by which nonlinearities in the equation of state work 
to break down the stratification of the Weddell Sea near Maud Rise? What are 
the relative roles of thermobaricity, cabbeling, and other NES issues? We will 
address these questions through detailed measurements of the small scale mixing 
processes in the upper ocean near Maud Rise and high resolution modeling of 
these processes. The goal is to produce parameterizations suitable for larger scale 


numerical models, and 


2.) How do the ocean circulation, ice and atmospheric forcing act to modulate and 
localize these processes? We will address this question by measuring the regional 
circulation and its variability near Maud Rise and by regional modeling of this 


circulation which incorporates the mixing parameterizations". 


During the MaudNESS deployment aboard the Research Vessel Ice Breaker 
(RVIB) Nathaniel B. Palmer, Antarctic Mesoscale Prediction System (AMPS) NWP 


2 


model forecast data were used in a simple ocean model to attempt to locate areas of 
possible deep convection. For the first time, this ocean model output was used as a 
planning tool for ship location and experimental purposes. Unfortunately, no deep 
convection was encountered during the MaudNESS deployment, as the Antarctic winter 


in 2005 was much milder than normal. 


This study was completed using archived AMPS data, ship observation data, and 
various seawater temperature and salinity profiles obtained during MaudNESS. Section 
If compares AMPS model output to observed conditions, both through time series and 
statistical analysis. Section III employs AMPS forecast, observed conditions, and 


variations of both as input to a simple one dimensional ocean model. 
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Figure 1. Antarctic region with Maud Rise’ area highlighted (After: 
https://www.cia.gov/cia/publications/factbook/reference_maps/antarctic.html). 


Hi. ATMOSPHERIC MODEL VERSUS OBSERVATIONS 


A. AMPS PMMS5 NUMERICAL WEATHER PREDICTION MODEL 


The numerical weather prediction (NWP) forecast model used during the 
MaudNESS deployment was the Antarctic Mesoscale Prediction System (AMPS). 
AMPS utilizes a polar version of the Pennsylvania State University-National Center for 
Atmospheric Research (NCAR) fifth generation mesoscale model (MM5). The MM5 
was modified for polar environments by the Polar Meteorology Group of the Byrd Polar 
Research Center at Ohio State University. Significant modifications made include 
improvements to cloud-radiation interaction, modified explicit ice phase microphysics, 
optimized boundary layer turbulence parameterization, sea ice surface type, and 


improved treatment of heat transfer through snow and ice surfaces. 


AMPS is initialized twice daily (00Z and 12Z) at the Mesoscale and Microscale 
Meteorology (MMM) Division of NCAR. The AMPS model is currently run in six 
domains (Figure 2). AMPS model data used on the MaudNESS deployment was from 
domain two, which at the time was run at a slightly coarser resolution of 30 kilometers. 
The model employs a terrain following staggered vertical grid, with thirty-two full sigma 
levels and thirty-one half sigma levels (half sigma levels are defined half way between 
full sigma levels). The lowest half sigma level is approximately thirteen meters above 
the surface (value used for bulk method used in this paper). The initial boundary 
conditions are driven by the National Center for Environmental Prediction's (NCEP) 
Global Forecasting System (GFS). Bromwich et al. (2005) state that: 

The observations available for assimilation include reports from 

radiosondes, surface SYNOP reports, automatic weather station (AWS) 

observations, ship reports, and buoys. Satellite-derived cloud-track winds 

are also assimilated in the 90-km grid. The system ingests sea ice data 


daily from the National Snow and Ice Data Center for its fractional sea ice 
depiction. 


Radiosonde data obtained during MaudNESS were received by MMM for 
assimilation into AMPS. It is unknown, however, if any data were missed or discarded 


through error checking procedures (Manning, 2007, personal communication). 


In May 2004 Three Dimensional Variational Data Assimilation (3DVAR) became 
the default analysis technique for AMPS MMS initialization. 


Domain \1 


Domain 3 
6.67 km 
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Antarctic Peninsula 





Figure 2. AMPS domains and horizontal resolutions. (After: 
http://www.mmm.ucar.edu/rt/mm5/amps/, 2007) 


B. OBSERVATIONAL WEATHER DATA 


The following shipboard meteorological data was collected every second, and 


averaged over ten second intervals: 
-Port and starboard wind speed (maximum, minimum, and average) 
-Port and starboard wind direction and standard deviation 
-Temperature, degrees Celsius, maximum and minimum 
-Relative humidity, single sample point, maximum and minimum 


-Barometric pressure (a five millibar bias was added to all observed surface 
pressure values due to instrument error. This error was verified by rawinsonde 


data). 
-Precision Spectral Pyranometer (PSP) (short wave radiation) 
-Precision Infrared Radiometer (PIR) (long wave radiation 


In addition, long wave and short wave radiation was measured using Eppley PIR 
and PSP sensors from the aft center of the helicopter deck (more accurate than ship PSP 
and PIR sensors, as they were cleaned and cleared of salt and ice daily). The sensors were 
mounted on gimbels so that they were level to approximately two degrees at all times. 
Data were collected every second and averaged into one minute records. Rawinsondes 


were also launched at a minimum of twice daily to collect upper air atmospheric profiles. 


C. COMPARISONS AND STATISTICS BETWEEN MODELED AND 
OBSERVED WEATHER DATA 


In order to achieve the stated purpose of this research, it was first necessary to 
determine how well the AMPS model performed. This was achieved by comparing 
archived AMPS model data with observational data gathered during the MaudNESS 
deployment. To ensure complete overlap between AMPS and observational data, the 
time period of comparison was limited to 1200Z 30 July 2005 to 1200Z 10 September 
2005 (Julian dates 211.5 to 253.5). During this time period the Research Vessel entered 


the ice on 30 July 2005, began ice egress on 4 September 2005 (Julian day 247) and 
exited the ice on 8 September 2005 (Julian day 251). 


The 12Z model run was sent to the Research Vessel daily and included forecasts 
in three hour time steps out to 120 hours. For this paper, continuous AMPS time series 
were created by combining daily forecast data beginning at 1200Z (time zero) and ending 
at O900Z the following day (time twenty-one). Ship meteorological data were collected 
every second and averaged over a ten second interval. For the chosen period of 
comparison this resulted in approximately 360 AMPS data points and almost 65,000 
observational data points. Using a simple MATLAB routine, observational data were 
merged with AMPS data (by Julian date) to provide a direct comparison between the two. 
Parameters chosen for comparison include surface pressure, near surface air temperature, 
relative humidity, specific humidity, wind speed and direction, downwelling long wave 


radiation, downwelling short wave radiation, latent heat flux, and sensible heat flux. 


Time series plots were then created for meteorological parameters, radiation, and 
heat fluxes for the entire period and weekly. Scatter plots were created comparing 


observed versus AMPS data for correlation analysis. 


For each parameter, a differential value was obtained by subtracting observation 
data values from AMPS model data values, and plotted as a distribution density. 
Statistical parameters such as minimum, maximum, mean, median, and root mean 
squared (RMS) error were calculated for AMPS data, observation data, and the 
differential data. In addition, correlation coefficient (R) and correlation coefficient 


squared (R’) were calculated for the differential data. 
1. Surface Pressure 


AMPS surface pressure was found to agree with observational relatively well. 
Time series plots (Figures 3 and 4) show AMPS over forecasting surface pressure on 
Julian days 220 to 223 and under forecasting surface pressure after Julian day 250. This 
can also be clearly seen on a scatter plot of observed versus AMPS surface pressure 


(Figure 5) and a differential distribution plot (Figure 6). 


Basic statistics for AMPS and observation surface pressure are shown in Table 1. 


Statistics run on differential values (AMPS minus observation) are shown in Table 2. 


Evaluation of statistical and correlation data indicate that AMPS was handling the 


synoptic pattern relatively well. 
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Surface pressure time series 30 July 2005 to 10 September 2005. Black line is 
observational data and red line is AMPS data. 
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equal to one drawn. Values for bias, differential root mean squared error, and 
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12 


Difference of Modeled and Observed 
Surface Pressure Distribution, 30Jul05 to 10Sep06 
90 T T T T T 





Number of Occurances 











a 
-30 -20 -10 Qa 10 20 30 
Surface Pressure Difference (mb) 


Figure 6. Distribution plot of AMPS surface pressure minus observation surface 


pressure. 
Nomen eedate Max | Mean | Median 
points sane 


AMPS Surface Pressure (mb) 953.2 | 1010.5 | 984.8 | 985.1 10.1 
Observation Surface Pressure (mb) 5 950.5 | 1017.6 | 985.6 | 987.5 12.6 


Table 1. Minimum, maximum, mean, median, and root mean squared values for AMPS 








derived surface pressure (row one) and observation surface pressure (row two). 
Number of Min Max MS Diff R R 
data points Error 
AMPS Surf Press minus 
Observation Surf Press 360 -28.6339 | 24.2523 | -0.8471 7.3879 0.656 | 0.4304 
(mb) 


Table 2. Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS minus observation surface 
pressure. 
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2 Near Surface Air Temperature 


As can be seen in Figures 7 and 8, near surface air temperature was quite often 
over forecast by AMPS. This can also be visualized in Figure 9 where a majority of the 
points fall above the correlation line, and in Figure 10 where the distribution density is 


skewed to the positive differential side. 


Basic statistics for AMPS and observation near surface air temperature are given 


in Table 3 and for differential data in Table 4. 


Air Temp Comparison 30Jul05 to 10Sep05 
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Figure 7. Near surface air temperature time series 30 July 2005 to 10 September 2005. 
Black line is observational data and red line is AMPS data. 
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Air Temp Comparison 30Jui05 to 06Aug0S 






































Air Temp Comparison 07Aug05 to 14Aug05 
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Julian Day Observation. black AMPS Nowcast...red 


Julian Day Observation..black AMPS Nowcast...red 


Near surface air temperature weekly time series. Black line is observational 
data and red line is AMPS data. a.) 30Jul05 to O6Aug05, b.) O07AugO5 to 
14Aug05, c.) 16Aug05 to 23Aug05, d.) 24Aug05 to 30Aug05, e.) 31Aug05 to 
06Sep05, and f.) O7Sep05 to 10Sep05. 


Figure 8. 
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Observed vs Modeled Near Surface Temperature, 
30Jul05 to 10Sep05 


Bias: 3.4952 °C 
Diff RMS Error: 4.2165 °C 
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Figure 9. Scatter plot of observed versus modeled near surface air temperature with 


correlation line equal to one drawn. Values for bias, differential root mean 
squared error, and correlation coefficient squared (R’) are shown in the upper 
right hand corner. 


16 


Difference of Modeled and Observed 
Near Surface Temperature Distribution, 30Jul05 to 105ep05 


Number of Occurances 








0 5 
Near Surface Temperature Difference (C) 


Figure 10. Distribution plot of AMPS near surface air temperature minus observation 
near surface air temperature. 


pins anes 


AMPS Air Temperature (°C) -21.5971 -11.862 | -11.8007 | 4.5652 





Observation Air Temperature 
CC) 
Table 3. Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived near surface air temperature (row one) and observation near surface air 
temperature (row two). 


-25.1 : -15.3572 
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Number of Min Max MS Diff R 

data points Error 
AMPS Air Temp minus 

Observation Air Temp (°C) -6.6809 | 14.5713 | 3.4952 | 4.2165 | 0.4469 | 0.1998 


Table 4. Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS minus observation near 
surface air temperature. 






2 Relative Humidity 


AMPS relative humidity was found to have no correlation to observed relative 
humidity, as can be seen in time series (Figures 11 and 12) and scatter plots (Figure 13). 
Based on this, differential distribution density (Figure 14) and basic statistical analysis 
(Tables 5 and 6), it was determined that AMPS had no skill in forecasting relative 


humidity. 


A comparison of relative humidity with respect to ice was also performed with 


similar results. Relative humidity with respect to ice was calculated as follows: 





RH,,, =RH| = (1) 


where RHice is relative humidity with respect to ice, RH is relative humidity, e, is 
saturation vapor pressure (both AMPS values and observed), and @gice is saturation vapor 


pressure with regards to ice (Both AMPS values and observed). 


Time series (Figure 15) showed that overall, observed values of relative humidity 
with respect to ice were at saturation up to Julian day 250. AMPS relative humidity with 
respect to ice did not follow this trend, with values as low as 65 % to values as high as 
125%. A scatter plot (Figure 16) of observed versus AMPS relative humidity with 
respect to ice showed no clear trends along a one-to-one correlation line. A differential 
distribution density plot (Figure 17) also showed wide variability between modeled and 
observed values, and a slight trend towards negative differences. Basic statistics and 


correlation are shown in Tables 7 and 8. Note that the correlation coefficient shown in 
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Table 8 is much less than | (0.0041) indicating no correlation between modeled and 


observed relative humidity with respect to ice. 


The final comparison was based on Andreas et al. (2002) which found that in a 
polar environment relative humidity with respect to ice was almost always near one 
hundred percent. With this in mind, equation (1) was rearranged to solve for RH with 
RHice set at one hundred percent (using AMPS derived saturation vapor pressure values). 
Time series comparisons of observed relative humidity and relative humidity calculated 
from relative humidity with respect to ice at saturation are shown in Figure 18. Relative 
humidity calculated from relative humidity with respect to ice followed observed relative 
humidity close than AMPS forecast relative humidity, but missed all of the drying trends 
seen in the observed values. The scatter plot shown in Figure 19 shows a slight increase 
in grouping around a one-to-one correlation line, and differential distribution density 
(Figure 20) shifted slightly away from negative differences. Basic statistics and 
correlation are shown in Tables 9 and 10. Correlation coefficient and correlation 
coefficient squared (0.2842 and 0.0808 respectively) indicate slightly better correlation 


when relative humidity with respect to ice is set to saturation than straight AMPS output. 
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RH Comparison 30Jul05 to 10Sep0S 
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Figure 11. Relative humidity time series 30 July 2005 to 10 September 2005. Black line 


is observational data and red line is AMPS data. 
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RH Comparison 30Jul0S to 06Aug0S 
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RH Comparison 31Aug05 to O6Sep0S 
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Julian Day Observation...black AMPS Nowcast...red 
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RH Comparison 07Aug05 to 14Aug05 
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RH Comparison 24Aug05 to 30Aug05 
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RH Comparison 07Sep05 to 10Sep0S 
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254 


Relative humidity weekly time series. Black line is observational data and red 
line is AMPS data. a.) 30Jul0S5 to 06Aug05, b.) O7Aug05 to 14Aug05, c.) 
16Aug05 to 23Aug05, d.) 24Aug05 to 30Aug05, e.) 31Aug05 to 06Sep05, and f.) 
07Sep05 to 10Sep05. 


Figure 12. 


21 





Bias: -0.2437 % 
Diff RMS Error: 8.3967 % 
R?: 0.0018 


Observed vs Modeled Relative Humidity, 
30Jul05 to 10Sep05 





110 - 
+ 
* a 
* a 
ee 
en 
100+ . % - - 7 
* * = 
+ a 
* +* * * a 
7 
* + 
* ** * 
* + 4+ . Pa 
+ ~ 
hs ». * : = et * 
+ + + *+* ee + 
* = ++ + + +<+ . * 
‘“ + * . * t" % * 
L * a 
90 * * o** a0 * * 4 ee : ye * + * 
+ *, At 
Ps ** *{ #E* tet + + * 4 * + 
& * * Os ert + + 
= * * e ~ +n * % *ee * * 
2 + ° * ate E* Ua ate OF ~ * 
3 * + tt * % * 
€ * + + ‘+ ~~ + #*¢ * 
5 Ps ee * 0° » © AE tte #9 ++ * * 
= ; + + eee + £ wtaee FFL, ‘ 
2 80+ w- # + Fee 
= * + KK «2 * #* +* 
o * a * «+ » * * & + 
a Pal » *, + bl 
3 *, . + * * ** 
2 SF * 
oF 7 + + 
3 a * ~  * Hy 
a 
= = 7 
a + 
70}- ge * + 
Pe + 
- bs 
we * 
y ‘ * 
+ 
. * 
60/- 
+ 
+ 
50 ! L L ! ! l } 
65 70 75 80 85 90 95 100 105 


Figure 13. 


Observed Relative Humidity (%) 


Scatter plot of observed versus modeled relative humidity with correlation line 
equal to one drawn. Values for bias, differential root mean squared error, and 
correlation coefficient squared (R’) are shown in the upper right hand corner. 
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Difference of Modeled and Observed 
Relative Humidity Distribution, 30Jul05 to 10Sep05 


Number of Occurances 











a 
Relative Humidity Difference (%) 


Figure 14. Distribution plot of AMPS relative humidity minus observation relative 
humidity. 


data points 


AMPS Relative Humidity (%) -21.5971 85.3236 | 85.2671 


Observation Relative Humidity 
(%) 





-25.1 85.5672 





Table 5. Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived relative humidity (row one) and observation relative humidity (row two). 
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Table 6. 


Relative Humidity wrt Ice (%) 


Figure 15. 


130 


120} 


90} 


80} 


70} 


AMPS Air Temp minus 
Observation Air Temp (°C) 










MS 
data points 
Error 
-28.0775 | 22.371 | -0.2437 | 8.3967 | 0.0427 | 0.0018 


Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS minus observation relative 
humidity. 


RH wrt Ice Comparison 30Jul05 to 10Sep0S 
i T T T T T T ] 
Observation 
—— AMPS 


PTR 





i 
215 220 225 230 235 240 245 250 255, 
Julian Day 


Relative humidity with respect to ice time series 30 July 2005 to 10 
September 2005. Black line is observational data and red line is AMPS data. 
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Bias: -3.4815 % 








Observed vs Modeled Relative Humidity wrt Ice, Diff RMS Error: 10.3727% 
30Jul0S to 10Sep05 
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Figure 16. Scatter plot of observed versus modeled relative humidity with respect to ice 


with correlation line equal to one drawn. Values for bias, differential root mean 
squared error, and correlation coefficient squared (R’) are shown in the upper 
right hand corner. 
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Difference of Modeled and Observed 
Relative Humidity wrt Ice Distribution, 30Jul05 to 10Sep05 
—— — 
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Relative Humidity Difference (%) 
Figure 17. Distribution plot of AMPS relative humidity with respect to ice minus 


observation relative humidity with respect to ice. 






uuniiyer . Min Max Mean Median RMS Error 
data points 


AMPS Relative 
Humidity of Ice(%) 65.5469 | 123.5048] 96.6991 95.8189 









Observation Relative 
Humidity of Ice(%) 75.378 | 107.5769} 100.1805 | 102.2276 
Table 7. Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived relative humidity with respect to ice (row one) and observation relative 


humidity with respect to ice (row two). 
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Table 


Relative Humidity (%) 


AMPS RH of Ice minus 
Ob RH of Ice (%) -39.3284 | 25.3316 |-3.4815] 10.3727 | 0.0041 | 1.72E-05 










Number of Min Max MS Diff R? 
data points Error 


8. Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS minus observation relative 
humidity with respect to ice. 


RH Comparison 30Jul05 to 10Sep05 
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Figure 18. Relative humidity time series 30 July 2005 to 10 September 2005. Black line 


is observation, red line is relative humidity calculated from relative humidity of 
ice at saturation. 
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Bias: 2.8049 % 








Observed RH vs RH wrt RH of Ice at Saturation, Diff RMS Error: 5.1327 % 
30Jul0S to 10Sep05 
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Figure 19. Scatter plot of observed relative humidity versus relative humidity calculated 


from relative humidity of ice at saturation, with correlation line equal to one 
drawn. Values for bias, differential root mean squared error, and correlation 
coefficient squared (R°) are shown in the upper right hand corner. 
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Difference of RH wrt RH of Ice at Saturation and Observed 


Relative Humidity Distribution, 30Jul05 to 10Sep05 
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Figure 20. Distribution plot of relative humidity calculated from relative humidity of ice 
at saturation minus observation relative humidity. 


Number of 
data points 


Relative Humidity 
wrt Relative 
Humidity of Ice at 
Saturation (%) 


80.3716 


99.1549 


88.3721 


88.3215 


RMS Error 





Observation Relative 
Humidity (%) 





67.03 








85.5672 





85.67 








Table 9. Minimum, maximum, mean, median, and root mean squared values for relative 
humidity calculated from relative humidity of ice at saturation (row one) and 


observation relative humidity (row two). 
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Number of { MS Diff 2 
[amon | Min | Max | ios |e | ® | 
RH wrt RH of Ice at 

Saturation minus 360 -6.2233 | 19.4923 | 2.8049 | 5.1327 | 0.2842 | 0.0808 
Observation RH (%) 


Table 10. | Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for relative humidity calculated from 
relative humidity of ice at saturation minus observation relative humidity. 






4. Specific Humidity 


Specific humidity, with a dependence on air temperature, followed a similar trend 
of being over forecast by AMPS. This can be visualized in time series (Figures 21 and 
22), scatter (Figure 23) differential distribution density (Figure 24) plots, and basic 
statistics (Tables 11 and 12). 
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Figure 21. Specific humidity time series 30 July 2005 to 10 September 2005. Black line 
is observational data and red line is AMPS data. 
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Specific Humidity Comparison 30.005 to O8Aug0S Specific Huridity Comparison 07 Aug0S to 14Aug05 
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Figure 22. Specific humidity weekly time series. Black line is observational data and red 


line is AMPS data. a.) 30Jul05 to 06Aug05, b.) O7AugO05 to 14Aug05, c.) 
16Aug05 to 23Aug05, d.) 24Aug05 to 30Aug05, e.) 31Aug05 to 06Sep05, and f.) 
07Sep05 to 10Sep0S5. 
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Figure 23. 





Observed vs Modeled Specific Humidity, 
30Jul0S to 10Sep05 
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Scatter plot of observed versus modeled specific humidity with correlation 
line equal to one drawn. Values for bias, differential root mean squared error, and 
correlation coefficient squared (R’) are shown in the upper right hand corner. 
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Difference of Modeled and Observed 
Specific Humidity Distribution, 30Jul05 to 10Sep05 
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Figure 24. Distribution plot of AMPS specific humidity minus observation specific 
humidity. 
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Table 11. | Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived specific humidity (row one) and observation specific humidity (row two). 


a2 






points Error 
AMPS Spec Humidity minus 
Ob Spec Humidity (g/kg) -1.1812} 1.9932 | 0.2899 | 0.4942 | 0.5321 | 0.2831 


Table 12. Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS minus observation 
specific humidity. 





5; Near Surface Wind Speed 


Near surface wind speed time series (Figures 25 and 26) showed that AMPS 
forecast values tended to agree well with observed values. The major error was found to 
be under-forecasting of wind speed (Julian days 215, 222 to 223, and 243 to 245). Wind 
speed under-forecast periods tend to occur when AMPS over-forecast surface pressure. 
A scatter plot of observed versus AMPS wind speeds (Figure 27), showed relatively good 
grouping of data points around the one-to-one correlation line. Differential distribution 
density is shown in Figure 28. Basic statistics and correlation are presented in Tables 13 


and 14. 


AMPS forecast wind speed was found to correlate relatively well (correlation 
coefficient of 0.551 shown in Table 14) with observed wind speed. This reinforced the 
idea that AMPS was handling the synoptic situation well, but slightly under-forecasting 


storm intensity. 
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Wind Speed Comparison 30Jul05 to 10Sep05 
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Figure 25. Near surface wind speed time series 30 July 2005 to 10 September 2005. 
Black line is observational data and red line is AMPS data. 
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Wind Speed Comparison 30.1005 to O6Aug05 Wind Speed Companson 07 Aug to 14Aug05 
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Figure 26. Near surface wind speed weekly time series. Black line is observational data 
and red line is AMPS data. a.) 30Jul05 to O6Aug05, b.) O7Aug05 to 14Aug05, c.) 
16Aug05 to 23Aug05, d.) 24Aug05 to 30Aug05, e.) 31 Aug05 to 06Sep05, and f.) 
07Sep05 to 10Sep05. 
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Figure 27. Scatter plot of observed versus modeled near surface wind speed with 


correlation line equal to one drawn. Values for bias, differential root mean 
squared error, and correlation coefficient squared (R’) are shown in the upper 
right hand corner. 
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Figure 28. Distribution plot of AMPS near surface wind speed minus observation near 
surface wind speed. 


umber of 
sry | Mm_| ton | nem [tote | oe 


AMPS Wind Speed (m/s) 20.4145 
Observation Wind Speed (m/s) . 29.35 
Table 13. | Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived near surface wind speed (row one) and observation near surface wind 


speed (row two). 
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Number of Min May MS Diff; R 
data points Error 
AMPS Wind Speed minus Ob 

Wind Speed (m/s) -11.7736 | 8.8707 |-0.7375 | 2.9507 | 0.551 | 0.3036 


Table 14. | Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS minus observation near 
surface wind speed. 


6. Near Surface Wind Direction 


Weekly time series plots (Figure 29, presented for better clarity) showed AMPS 
near surface wind direction compared well with observed wind direction, with notable 


exceptions on Julian days 216, 219 to 220, 242 to 244, and 252. 


Figure 30 presents a scatter plot of observed versus AMPS wind direction, and 
Figure 31 shows differential distribution density. It should be noted that these plots and 
basic statistics (Table 15) can be misleading, appearing to have large error when wind 
directions are northerly. This is due to the wind direction being represented from 0 to 
360 degrees clockwise, with north being at the 0 and 360 degree point. With this in 
mind, and inspection of the time series plots shows that wind direction was indeed 


forecast well by AMPS. 
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Near surface wind direction weekly time series. Black line is observational 
data and red line is AMPS data. a.) 30Jul05 to O6Aug05, b.) O7AugO5 to 
14Aug05, c.) 16Aug05 to 23Aug05, d.) 24Aug05 to 30Aug05, e.) 31Aug05 to 
06Sep05, and f.) O7Sep05 to 10Sep05. 


Figure 29. 
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Observed vs Modeled Near Surface Wind Direction, 
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Figure 30. 


Observed Near Surface Wind Direction (degrees) 


Scatter plot of observed versus modeled near surface 


correlation line equal to one drawn. 
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wind direction with 


Difference of Modeled and Observed 
Near Surface Wind Direction Distribution, 30Jul05 to 10Sep05 
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Figure 31. Distribution plot of AMPS near surface wind direction minus observation near 
surface wind direction. 


Number of data points 


AMPS Wind Direction (degrees) 215.1103 | 229.5342 80.01 
Observation Wind Direction (degrees) 214.8948 : 81.6885 
Table 15. | Mean, median, and root mean squared values for AMPS derived near surface 
wind direction (row one) and observation near surface wind direction (row two). 








7, Downwelling Long Wave Radiation 


AMPS long wave radiation forecasts were found to have no correlation with 


observed. Time series plots are shown in Figures 32 and 33, scatter plot in Figure 34, and 
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differential distribution density plot in Figure 35. As can be seen in Figure 34 AMPS had 
no skill in forecasting long wave radiation. Basic statistics and correlation calculations 


(Tables 16 and 17) confirm this. 


Longwave Radiation Comparison 30Jul05 to 10Sep05 
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Figure 32. Downwelling long wave radiation time series 30 July 2005 to 10 September 


2005. Black line is observational data and red line is AMPS data. 
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Longwave Radiation Comparison 30Jui0S to O6Aug05 
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Longwave Rackation Comparison 31Aug05 to O6Sep05 











Longwave Radiation Comparison 07Aug0S to 14Aug05 
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Longwave Radiation Comparison 24Aug05 to 1OAug05S 
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Longwave Radation Comparison 07Sep05 to 10Sep05 
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Downwelling long wave radiation weekly time series. Black line is 
observational data and red line is AMPS data. a.) 30Jul05 to 06Aug0S5, b.) 
07Aug05 to 14Aug05, c.) 1l6AugO05 to 23Aug05, d.) 24Aug05 to 30Aug05, e.) 
31Aug05 to 06Sep05, and f.) O7Sep05 to 10Sep05. 





Bias: -13.8245 W/m? 
Diff RMS Error: 51.4731 W/m? 


Observed vs Modeled Downwelling Longwave Radiation, 
30Jul05S to 10Sep05 
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Figure 34. Scatter plot of observed versus modeled downwelling long wave radiation 


with correlation line equal to one drawn. Values for bias, differential root mean 
squared error, and correlation coefficient squared (R’) are shown in the upper 
right hand corner. 


45 


Difference of Modeled and Observed 
Downwelling Longwave Radiation Distribution, 30Jul05 to 10Sep05 
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Figure 35. Distribution plot of AMPS downwelling long wave radiation minus 


observation downwelling long wave radiation. 






RMS 
Error 






Number of | Min Max Mean | Median 
data points 


AMPS Downwelling Long Wave 
Radiation (W/m?) 106.0907 | 296.1095 | 201.7034 | 205.1429 | 41.9733 









Observation Downwelling Long 
Wave Radiation (W/m?) 136.6 316.6 | 215.5279} 223.5 45.6767 


Table 16. Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived downwelling long wave radiation (row one) and observation downwelling 


long wave radiation (row two). 
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points Error 
AMPS LW Radn minus 
Ob LW Radn (W/m?) - 122.389 | 120.0858 | -13.8245 | 51.4731 | 0.0977 | 0.0095 


Table 17. | Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS minus observation 
downwelling long wave radiation. 





8. Downwelling Short Wave Radiation 


Figures 36 and 37 show good overlap of AMPS forecast downwelling short wave 
radiation compared to observed downwelling short wave radiation. One feature to note in 
these figures are several times where AMPS forecast negative values for short wave 
radiation (particularly Julian days 258 to 251). It is unknown why this occurred and was 
not corrected for in statistical calculations (Tables 18 and 19). A scatter plot of observed 
versus AMPS forecast downwelling radiation is shown in Figure 38. Differential 
distribution density (Figure 39) shows the majority of differences are slightly positive. 
These plots, as well as a correlation coefficient of 0.7891 (Table 19), tend to indicate 


AMPS showed high skill in forecasting downwelling short wave radiation. 
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Shortwave Radiation Comparison 30Jul05 to 10Sep05 
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Figure 36. Downwelling short wave radiation time series 30 July 2005 to 10 September 
2005. Black line is observational data and red line is AMPS data. 
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SPorteave Radation Comparison 10Ju05 to 06Aug05 Shortwave Radation Comparison O7 Aug0S to 14Aug05 
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Figure 37. Downwelling short wave radiation weekly time series. Black line is 


observational data and red line is AMPS data. a.) 30Jul05 to 06Aug0S5, b.) 
07Aug05 to 14Aug05, c.) 16Aug05 to 23Aug05, d.) 24Aug05 to 30Aug05, e.) 
31Aug05 to 06Sep05, and f.) O7Sep05 to 10Sep05. 
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Observed vs Modeled Downwelling Shortwave Radiation, Diff RMS Error: 33.9058 W/m2 
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Figure 38. Scatter plot of observed versus modeled downwelling short wave radiation 


with correlation line equal to one drawn. Values for bias, differential root mean 
squared error, and correlation coefficient squared (R*) are shown in the upper 
right hand corner. 
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Difference of Modeled and Observed 
Downwelling Shortwave Radiation Distribution, 30Jul05 to 10Sep0S 
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Distribution plot of AMPS downwelling short wave radiation minus 
observation downwelling short wave radiation. 











Number of 


data points mn 


-26.5966 


378.1499 


Mean 


42.1284 


Median 


RMS Error 


73.4132 





Observation Downwelling 


Short Wave Radiation -6.963 

(W/m?) 
Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived downwelling short wave radiation (row one) and _ observation 
downwelling short wave radiation (row two). 





Table 18. 
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358.6 





33.1134 





68.8514 









Number of Min Max MS Diff R 
data points Error 
AMPS SW Radn minus 

Ob SW Radn (W/m?) - 133.353 | 248.7499 | 9.015 33.9058 | 0.7891 | 0.6228 


Table 19. Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS minus observation 
downwelling short wave radiation. 





9, Latent and Sensible Heat Fluxes 


The AMPS forecast model provides direct latent and sensible heat flux outputs. 
For this research observed heat fluxes were calculated using the bulk method. In 
addition, AMPS bulk parameters (relative humidity, wind speed, air temperature, and 
surface pressure) were used to calculate modeled heat fluxes by the bulk method. Time 
series plots (Figures 40, 41, 46, and 47) display AMPS direct heat flux output (red), 
observed bulk method (black) and AMPS bulk method (cyan). As can be seen in these 
plots, bulk method values tend to be similar, while model output appears to have no 
correlation at all. The problem with this evaluation is that for the bulk methods, surface 
temperature for ice is assumed to be air temperature (no direct surface temperature 


measurements were obtained). 


Ice and ocean heat fluxes calculated by bulk methods were then combined into a 
total flux using a weighted average based on percent ice cover and open ocean (85% and 


15% respectively). For example, the calculation for sensible heat flux was: 
SNF rat = (85 )(shf,.) + CIS)ShA, oan ) (2) 


where shf, shf,., and shf,., are total, ice, and ocean sensible heat fluxes 


otal ? 
respectively. When assuming surface temperature is equal to near surface air 
temperature, all of the sensible heat flux (and most of the latent) are derived from the 
open ocean. Another issue was the percentages of ice cover and open ocean were kept 
constant, when in reality they were variable. Direct model output of heat fluxes showed 


more variability than bulk methods because surface temperature was not fixed to near 
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surface air temperature. Bulk method calculations were conducted using a program 


written by Dr. Peter Guest of the Naval Postgraduate School Department of Meteorology. 


Scatter plots were created for observed bulk method versus direct model output 
(Figures 42 and 48) and observed bulk method versus model bulk method (Figures 43 
and 49). As expected, observation versus direct AMPS output indicated no skill, while 
observed versus model bulk methods showed a correlation. Differential distribution 
density plots (Figures 44, 45, 50, and 51) and basic statistics and correlations (Tables 20, 
21, 22 and 23) show the same results. 
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Latent Heat Flux Comparison 30Jul05 to 10Sep05 
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Figure 40. Latent Heat Flux time series 30 July 2005 to 10 September 2005. Black line is 
observational data, red line is AMPS latent heat flux output, and cyan line is latent 
heat flux calculated by bulk methods using AMPS model output for bulk 
parameters. 
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Figure 41. Latent heat flux weekly time series. Black line is observational data, red line 


is AMPS latent heat flux output, and cyan line is latent heat flux calculated by 
bulk methods using AMPS model output for bulk parameters. a.) 30Jul05 to 
06Aug05, b.) O7Aug05 to 13Aug05, c.) 14Aug05 to 20Aug05, d.) 21Aug05 to 
27Aug05, e.) 28Aug05 to 03Sep05, and f.) 04Sep05 to 10Sep05. 
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Observed vs Modeled (Direct AMPS Output) Latent Heat Flux, Diff RMS Error: 8.7511 Wim? 
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Figure 42. Scatter plot of observed versus modeled latent heat flux (direct AMPS output) 


with correlation line equal to one drawn. Values for bias, differential root mean 
squared error, and correlation coefficient squared (R’) are shown in the upper 
right hand corner. 
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Observed Latent Heat Flux (W/m?) 


Scatter plot of observed versus modeled latent heat flux (calculated by bulk 
method using AMPS model output for bulk parameters) with correlation line 
equal to one drawn. Values for bias, differential root mean squared error, and 
correlation coefficient squared (R’) are shown in the upper right hand corner. 
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Difference of Modeled (Direct AMPS Output) and Observed 
Latent Heat Flux Distribution, 30Jul05 to 10Sep05 
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Figure 44. Distribution plot of AMPS latent heat flux (direct model output) minus 
observation latent heat flux. 
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Difference of Modeled (Bulk Method) and Observed 
Latent Heat Flux Distribution, 30Jul05 to 10Sep05 
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Latent Heat Flux (Wim?) 
Figure 45. Distribution plot of AMPS latent heat flux (calculated by bulk method using 
AMPS model output for bulk parameters) minus observation latent heat flux. 


Number of 


data points RMS Error 


AMPS Latent Heat Flux, 
Model Output (W/m’) 
AMPS Latent Heat Flux, 
Bulk Method (W/m2) 
Observation Latent Heat 
Flux (W/m’) 
Table 20. Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived (direct model output) latent heat flux (row one), AMPS derived 
(calculated by bulk method using AMPS model output for bulk parameters) latent 
heat flux (row two) and observation latent heat flux (row three). 


-11.6854 | 37.4711 ) 13.5072 | 13.2429 10.8533 





1.2293 | 26.6983 | 12.3249 12.4332 4.841 





-6.7852 | 36.2324 | 14.1492 | 14.7403 6.5025 
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Number of i MS Diff 
data points Error 


AMPS LHF minus Ob 
LHF (W/m’) 


-20.6235 | 22.3135 





AMPS Bulk Method 
LHF minus Ob LHF -21.3666 | 12.5859 
(W/m2) 
Table 21. | Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS (direct model output) 
minus observation latent heat flux (row one), and AMPS (bulk method) minus 

observation latent heat flux (row two). 
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Figure 46. Sensible heat flux time series 30 July 2005 to 10 September 2005. Black line 


is observational data, red line is AMPS sensible heat flux output, and cyan line is 
sensible heat flux calculated by bulk methods using AMPS model output for bulk 
parameters. 
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Sensible Heat Faux Comparison 30.2005 to O6Aug05 
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Sensible heat flux weekly time series. Black line is observational data, red line 
is AMPS sensible heat flux output, and cyan line is sensible heat flux calculated 
by bulk methods using AMPS model output for bulk parameters. a.) 30Jul05 to 
06Aug05, b.) O7Aug05 to 13Aug05, c.) 14Aug05 to 20Aug05, d.) 21Aug05 to 
27Aug05, e.) 28Aug05 to 03Sep05, and f.) 04Sep05 to 10Sep05. 


Figure 47. 
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Observed vs Modeled (Direct AMPS Output) Sensible Heat Flux, 


Bias: 1.332 Wim? 
Diff RMS Error: 23.320 W/m2 
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Figure 48. Scatter plot of observed versus modeled sensible heat flux (direct AMPS 


output) with correlation line equal to one drawn. Values for bias, differential root 
mean squared error, and correlation coefficient squared (R’) are shown in the 
upper right hand corner. 
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Bias: -4.9932 W/m? 





Observed vs Modeled (Bulk Method) Sensible Heat Flux, Diff RMS Error: 10.0669 W/m? 
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Figure 49. Scatter plot of observed versus modeled sensible heat flux (calculated by bulk 


method using AMPS model output for bulk parameters) with correlation line 
equal to one drawn. Values for bias, differential root mean squared error, and 
correlation coefficient squared (R) are shown in the upper right hand corner. 
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Difference of Modeled (Direct AMPS Output) and Observed 
Sensible Heat Flux Distribution, 30Jul05 to 10Sep05 
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Figure 50. Distribution plot of AMPS sensible heat flux (direct model output) minus 
observation sensible heat flux. 
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Difference of Modeled (Bulk Method) and Observed 
Sensible Heat Flux Distribution, 30Jul05 to 10Sep05 
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Figure 51. Distribution plot of AMPS sensible heat flux (calculated by bulk method 
using AMPS model output for bulk parameters) minus observation sensible heat 
flux. 


Number of : Median 
data points 


AMPS Sensible Heat Flux, 
Model Output (W fm?) -48.2254 | 94.7855 | 26.3397 | 26.402 | 27.3168 


AMPS Sensible Heat Flux, Bulk 
Method (W/m2) -6.0896 | 47.6757 | 20.2441 | 21.1381 9.3373 


inca Tal Bere -11.6588 | 58.8124 | 25.0077 | 25.5618 | 13.5801 
Table 22. Minimum, maximum, mean, median, and root mean squared values for AMPS 
derived (direct model output) sensible heat flux (row one), AMPS derived 

(calculated by bulk method using AMPS model output for bulk parameters) 


sensible heat flux (row two) and observation sensible heat flux (row three). 
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Number of MS Diff 
data points Error 


AMPS SHF minus Ob 
SHF (W/m?) -77.929 | 61.3296 23.32 
AMPS Bulk Method SHF 
minus Ob SHF (W/m2) -50.0183 | 23.5586 10.0576 


Table 23. Minimum, maximum, bias, differential mean squared error, correlation coefficient 
(R), and correlation coefficient squared (R’) for AMPS (direct model output) 
minus observation sensible heat flux (row one), and AMPS (bulk method) minus 
observation senisble heat flux (row two). 





























D. COMPARISON CONCLUSIONS 


AMPS surface pressure and relative humidity forecasts were found to have no 
skill when compared to observed values. It is believed this was due to AMPS handling of 
cloud and moisture parameters. It was also found that a slightly better relative humidity 
forecast could be obtained by simply setting relative humidity with respect to ice to 


saturation. 


The skill of AMPS in predicting latent and sensible heat fluxes is unknown due to 
the broad assumptions made in calculating observed values with the bulk method. When 
AMPS bulk parameter inputs were used to calculate fluxes with the bulk method, the 
results were found to be similar with observed values (even with AMPS' lack of skill in 
forecasting relative humidity) because of the dependence of latent heat flux on specific 


humidity difference and wind speed. 


AMPS near surface air temperature, and therefore specific humidity, were found 
to be over-forecast part of the time. The remaining parameters (surface pressure, wind 
speed and direction, and short wave radiation) performed well when compared to 
observed values. Surface pressure, wind speed and direction show that AMPS had a 


good performance in representing the synoptic situation. 
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When time series plots for near surface wind speed were compared with time 
series for surface pressure, it was noted that times when AMPS under-forecast wind 
speed tended to correspond to times when it over-forecast surface pressure. This led to 


the belief that AMPS also tended to slightly under-forecast the intensity of storms. 


Residual scatter plot analysis was also performed for all parameters. This was 
accomplished by taking the difference of the value of a modeled parameter minus the 
observed value and comparing this difference to all other meteorological parameters 
(AMPS and observed) on a scatter plot. No significant correlations were found during 


this process. 
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Ht. OCEAN MODEL 


A. 1-D OCEAN MODEL 


The ocean model used for this research is a simple one dimensional, ocean model 
developed by Dr. Miles McPhee. McPhee (2000) gives a detailed description of the 
model, which will not be restated here. Basically, the model ocean structure is initialized 
by a temperature-salinity profile (conductivity-temperature-depth (CTD) profiler cast 
data used for this research). Time series were created for surface friction velocity, ice 
heat conduction, and downwelling short wave radiation. These fields were then used to 


force the initialized conditions over a selected period of time. 


It is important to note that only one CTD profile could be used to initialize the 
model for the entire run, regardless of the length of the model run. The actual ocean 
structure would not remain exactly the same over time, and, the research vessel did not 
remain in the location of the CTD profile cast. Other limitations of this model presented 
in McPhee (2000) are divergence of ice cover and Ekman transport in the mixed layer are 
neglected, and any assumption of horizontal homogeneity implicit in one dimension 


modeling, particularly in the Maud Rise region, are suspect. 


These limitations would severely degrade ocean modeling requiring a high degree 
of accuracy. For the purposes of this research, modeling the ocean structure over the 
Maud Rise was not the goal. Sensitivity to changes in atmospheric forcing was desired, 


and the one-dimensional ocean model proved sufficient in that task. 


B. COMPILATION OF DATA FOR USE IN MODEL 


1. Model Constants 


The following is a list of constants and their values used in the ocean model runs: 
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a. Model Depth 


A model depth of 250 meters was chosen (pycnocline depth observed in 


CTD data was between 50 to 150 meters). 
b. Undersurface Roughness Length (Z,) 


Undersurface roughness length of 1.0 mm was used (McPhee 2007, 


personal communication). 
Cc. Initial Ice Thickness 
Initial average ice thickness over the Maud Rise area was set to 40 cm. 
d. Heat transfer coefficient 


Based on McPhee et al. (1999), a dimensionless heat transfer coefficient 


(Stanton number based on friction speed instead of mean speed) of 5.6 x 10° was used. 
e. Latitude 


Latitudes for model runes were the latitude of the CTD profile used to 


initialize the run. 
f Buoyancy Frequency Squared (N’) Limit 


Based on Mcphee (2000), a buoyancy frequency squared (N°) limit of 
1.5x10° s” was used. The pycnocline depth is defined as the depth in the water column 


where this value is first exceeded. 
2. Model Variables 


The following is a list that describes the variables used in the ocean model and 


how they were compiled. 
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a. Temperature/Salinity Profiles 


Temperature and salinity profiles were obtained from three separate CTD 
downcast profiles at stations 1, 65, and 92 (reason for these cast selections is given 
below). Depth, primary temperature, and primary salinity data were taken from the casts 


and compiled into text files for input into the ocean model. 
b. Downwelling Solar Radiation 


Two downwelling solar radiation time series files were created from 
AMPS model output and observation data (AMPS downwelling short wave radiation and 


observed downwelling short wave radiation). 
c. Friction Speed (ux) 


McPhee et al. (1999) characterize friction speed as the square root of 
kinematic stress at the interface. Using this definition, friction speed was calculated 


using the following equation: 


Uu, = Total ( 3 ) 


where u., is friction speed in meters/second, 1,,,,, 18 wind stress in Newton/meters” (note 


that the wind stess is total wind stress over ice and seawater, taking into account the 


percent of ice cover), and p,, is the density of seawater (1025 kilograms/meter’ (kg/m’) 
used). 


Four separate wind stress time series were created using bulk methods for 


the following conditions: 
-AMPS forecast wind speed 
-Observational wind speed 
-AMPS wind speed increased by a factor of two, and 
-Observational wind speed increased by a factor of two. 


a. 


(Weather conditions during MaudNESS were found to be unusually mild 
for a winter season in the Weddell Sea. For this reason, AMPS and observed wind 
speeds were increased by a factor of two in order to simulate more typical winter forcing 


conditions). 


d. Ice Conductive Heat Flux ( 9 ) 


Ice conductive heat flux is defined by McPhee (1992) as heat conduction 
in ice divided by the product of density and specific heat of seawater. Using this 


definition, ice conductive heat flux was calculated using the following equation: 


a (4) 
(a )(¢.)} 
where q is ice conductive heat flux (meters-Celsius/second), q,., is heat conduction in 
ice (with respect to ice temperature) in Watts/meter’, P,, 18 the density of seawater (1025 
kg/m? used), and c »,, 18 the specific heat capacity of seawater (Joules/kilogram-Celsius) 
at the seawater-ice interface. 


Bulk methods were used to create heat conduction in ice time series for 
the same four scenarios described above. Specific heat capacity of seawater at the 
seawater-ice interface was calculated using interface values from the appropriate CTD 


profiles. 
C. CTD STATIONS CHOSEN FOR MODEL RUNS 


Figure 52 shows the bathymetry of the Maud Rise area and all the locations of the 
MaudNESS deployment Phase I CTD stations. The three highlighted CTD stations were 


the locations chosen for this study. 


AZ 


The temperature-salinity profiles were used to initialize the ocean model at the 
latitude and Julian day of the CTD cast (summarized in Table 24). The model was then 
run for seven days from the date of initialization, for each of the scenarios listed in the 


previous section. 
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Figure 52. Maud Rise bathymetry and MaudNESS CTD stations. Location of CTD's 


used for ocean model initialization annotated. (After: Figure 8, MaudNESS Cruise 
Report) 


CTD Station Number 


Date 


Time (Z) 


Lat (dec) 


Long (dec) 





1 


2-Aug-05 


12:43:14 


-63.236 


0.01 





65 


9-Aug-05 


5:09:28 


-64.8 


2.88 





92 





Table 24. 





21-Aug-05 





16:07:37 


73 





-65.515 





1.312 


Date, time, latitude, and longitude data for CTD stations used in this study. 


1. CTD Station 1 


CTD station 1 was selected as the station farthest away from the Maud Rise. The 
water column at this location was relatively stable when compared to the other two 
locations. The potential density difference across the pycnocline was approximately 
0.175 kg/m’, as can be seen in Figure 53. There were also relatively large differences in 
potential temperature and salinity across the pycnocline (approximately 2.8 degrees 


Celsius (°C) and 0.375 practical salinity units (psu) respectively). 


CTD Station 1, Downcast 
Potential Ternperature Salinity Potential Density 
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Figure 53. Initial potential temperature, salinity, and potential density profiles for CTD 


station 1 downcast. 


De CTD Station 65 


CTD station 65 was chosen due to the location over the Maud Rise in the Taylor 


column. This water column exhibited marginal stability with a potential density jump of 
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0.08 kg/m? (illustrated in Figure 54). The potential temperature and salinity differences 
are also smaller than at station 1 (approximately 2 °C and 0.2 psu respectively). 


CTD Station 65, Downcast 
Potential Temperature Salinity Potential Density 
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Figure 54. Initial potential temperature, salinity, and potential density profiles for CTD 


station 65 downcast. 


3. CTD Station 92 


CTD station 92 was located to the southwest of the Maud Rise in an area that 
Muench et al. (2001) called the warm pool. They further define this area as a 
"dynamically necessary region of positive (cyclonic) vorticity that is associated with a 
Taylor column over the rise". It was the least stable of the three locations chosen for this 
study. The potential density difference (Figure 55) at this location was approximately 
0.05 kg/m*. The potential temperature difference was approximately 2.3 °C (slightly 
higher than at station 65 but expected in the warm pool), and the salinity difference was 


approximately 0.18 psu. 
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CTD Station 92, Downcast 
Potential Temperature Salinity Potential Density 
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Figure 55. Initial potential temperature, salinity, and potential density profiles for CTD 


station 92 downcast. 


D. COMPARISON OF OCEAN MODEL RESULTS FORCED BY AMPS AND 
OBSERVATIONAL DATA 


1. CTD Station 1 


Figure 56 shows friction speed for the four different scenarios run. With the 
exception of day four, observed forcing was shown to result in a slightly higher friction 
speed than AMPS forcing. This result carried over to the friction speeds when both 


AMPS and observed wind speeds were increased. 


Figure 57 shows temperature fluxes for all four scenarios. Both cases of 
increased wind speed (Figure 57 c. and d.) indicate relatively strong and deep positive 


temperature flux when compared to AMPS forecast and observed values (Figure 57 a. 
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and b.). For AMPS forecast and observation inputs, a relatively large negative salinity 
flux persisted over the entire period (Figure 58 a. and b.). For increased wind speed 
scenarios, a large positive salinity flux was initially present with small intrusions of 
negative salinity flux later in the period (Figure 58 c. and d.). Buoyancy flux (Figure 59) 


was found to closely follow salinity flux. 


For the AMPS forecast and observation input, the weaker temperature flux 
allowed the sea ice to grow, which in turn caused negative salinity flux (increased salinity 
in the upper layer due to brine rejection of the freezing seawater). Figure 60 shows ocean 


model predicted ice thickness (initially 40 cm) increasing almost linearly after day one. 


When forcing was increased, the ocean model predicted an ice thickness decrease 
(Figure 61), with small increases in thickness around day three and days six and seven 
(corresponding to negative buoyancy and salinity fluxes at that time). The increased 
temperature flux for these runs immediately caused the sea ice to melt, thereby mixing 
fresh water with the upper layer water and creating positive salinity and buoyancy fluxes 


in the upper layer. 


Potential temperature, salinity, and potential density profiles are shown in Figures 
62, 63, and 64 respectively. Shown are initial CTD profiles (blue), and ocean model 
prediction at the conclusion of the run (with AMPS input (black), observed input (red), 
AMPS wind speed increased by factor of two (green), and observed wind speed increased 
by factor of two (cyan)). These profiles agree with results obtained from contour and ice 


thickness plots discussed above. 


For CTD station 1, the ocean model forced by observed conditions led to the least 
stable profile after seven days. Under this scenario, potential density difference across 
the pycnocline decreased to approximately 0.10 kg/m*, compared to an initial value of 


0.175kg/m°. 
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Figure 56. Friction speed (ux) contour plots from surface to 250 meters at CTD station 1 
for Julian days 214 to 221: a.) AMPS forecast, b.) Observed conditions, c.) AMPS 
forecast wind speed increased by factor of two, and d.) Observed wind speed 
increased by factor of two. 
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Figure 57. 
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Temperature flux (<w'T'>) contour plots from surface to 250 meters at CTD 


station | for Julian days 214 to 221: a.) AMPS forecast, b.) Observed conditions, 
c.) AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two. 
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Figure 58. Salinity flux (<w'S'>) contour plots from surface to 250 meters at CTD station 


1 for Julian days 214 to 221: a.) AMPS forecast, b.) Observed conditions, c.) 
AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two (Note white shaded contours are areas of 


relatively large negative values). 
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Figure 59. 
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Buoyancy flux (<w'b'>) contour plots from surface to 250 meters at CTD 


station | for Julian days 214 to 221: a.) AMPS forecast, b.) Observed conditions, 
c.) AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two (Note white shaded contours are areas of 


relatively large negative values). 
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Figure 60. Ocean model predicted ice thickness at CTD station 1, Julian days 214 to 221, 
for a.) AMPS forecast, and b.) Observed wind forcing conditions. 
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Hee: CTD 1, AMPS Wind Speed Increased by Factor of 2 
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Figure 61. Ocean model predicted ice thickness at CTD station 1, Julian days 214 to 221, 
for a.) AMPS forecast wind speed increased by factor of two, and b.) Observed 
wind speed increased by factor of two. 
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Figure 62. Potential temperature profiles for CTD station 1. Profiles shown are CTD 


initial (blue), and ocean model prediction seven days later (with AMPS forcing 
(black), observed forcing (red), AMPS wind speed increased by factor of two 
(green), and observed wind speed increased by factor of two (cyan)). 
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Figure 63. Salinity profiles for CTD station 1. Profiles shown are CTD initial (blue), and 


ocean model prediction seven days later (with AMPS forcing (black), observed 
forcing (red), AMPS wind speed increased by factor of two (green), and observed 
wind speed increased by factor of two (cyan)). 
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Figure 64. Potential density profiles for CTD station 1. Profiles shown are CTD initial 
(blue), and ocean model prediction seven days later (with AMPS forcing (black), 
observed forcing (red), AMPS wind speed increased by factor of two (green), and 
observed wind speed increased by factor of two (cyan)). 


2. CTD Station 65 


Figure 65 shows friction speed for the four different scenarios run at CTD station 
65 (Taylor column over the Maud Rise). Friction speed values between AMPS forecast 
and observed conditions (Figure 65 a. and b.) were very similar (observed conditions 
slightly larger). Again, when AMPS and observed forcing were increased (Figure 65 c. 


and d.), friction speed increased and high values extended deeper into the water column. 


Relatively large positive temperature flux (note the temperature flux contour 
values are close to two orders of magnitude larger than CTD station 1) for the AMPS 
forecast scenario (Figure 66 a.) occurred late in the model run (end of days five and 


seven). For observed conditions (Figure 66 b.), positive temperature fluxes occur 
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between days two and three, three and four, five and six, and six and seven. For each of 
the increased forcing scenarios (Figure 66 c. and d.), a large positive temperature flux 


area appeared between days two and three. 


Salinity and buoyancy fluxes (Figures 67 and 68 respectively) tended to be 
negative the entire period for all scenarios. A notable exception to this is the large area of 
positive salinity and buoyancy fluxes (Figures 67 c. and d. and 68 c. and d.) from 
approximately 75 to 150 meters between days two and three for the increased forcing 
scenarios (this area of positive salinity and buoyancy fluxes corresponds to the large 


positive temperature flux discussed in the previous paragraph). 


Figures 69 and 70 show an almost linear increase in ocean model predicted ice 
thickness for all four scenarios (note in Figure 70 a. and b., the rate of ice growth 


decreases around day three in agreement with temperature, salinity, and buoyancy flux). 


Potential temperature, salinity, and potential density profiles are shown in Figures 
71, 72 and 73, respectively. Initial CTD profiles and ocean model prediction at the 
conclusion of the run are as described for CTD station 1. As with CTD station 1, profiles 


agree well with contour and ocean model predicted ice thickness plots. 


At CTD station 65, the increased observed forcing scenario resulted in the least 
stable water column at the end of the seven day model run. Potential density difference 
across the pycnocline (Figure 73) decreased from an initial value of 0.08 kg/m? to 


approximately 0.02 kg/m’. 


87 


u,, CTD 65, AMPS Input u,, CTD 65, Observation Input 


0 
| =~ 















€ 
2801 222 223 224 225 226 227 228 222 223 224 225 226 227 228 
Day of 2005 Day of 2005 
u,, CTD 65, AMPS Wind Speed Increased by Factor of 2 u,, CTD 65, Observation Wind Speed Increased by Factor of 2 
it) 0 
€ € 
-200 
201 222.0 «-223—C—224 225s 226 227) 228 nae-ey 222 42230 «6224 «= 225, 226 227) 228 
Day of 2005 Day of 2005 
Figure 65. Friction speed (u«) contour plots from surface to 250 meters at CTD station 65 


for Julian days 221 to 228: a.) AMPS forecast, b.) Observed conditions, c.) AMPS 
forecast wind speed increased by factor of two, and d.) Observed wind speed 
increased by factor of two. 
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Figure 66. Temperature flux (<w'T'>) contour plots from surface to 250 meters at CTD 
station 65 for Julian days 221 to 228: a.) AMPS forecast, b.) Observed conditions, 
c.) AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two. 
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Salinity flux (<w'S'>) contour plots from surface to 250 meters at CTD station 
65 for Julian days 221 to 228: a.) AMPS forecast, b.) Observed conditions, c.) 
AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two (Note white shaded contours are areas of 
relatively large negative values). 
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Figure 68. Buoyancy flux (<w'b'>) contour plots from surface to 250 meters at CTD 


station 65 for Julian days 221 to 228: a.) AMPS forecast, b.) Observed conditions, 
c.) AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two (Note white shaded contours are areas of 
relatively large negative values). 
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Figure 69. Ocean model predicted ice thickness at CTD station 65, Julian days 221 to 
228, for a.) AMPS forecast, and b.) Observed wind forcing conditions. 
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Figure 70. Ocean model predicted ice thickness at CTD station 65, Julian days 221 to 
228, for a.) AMPS forecast wind speed increased by factor of two, and b.) 
Observed wind speed increased by factor of two. 
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Figure 71. Potential temperature profiles for CTD station 65. Profiles shown are CTD 
initial (blue), and ocean model prediction seven days later (with AMPS forcing 
(black), observed forcing (red), AMPS wind speed increased by factor of two 
(green), and observed wind speed increased by factor of two (cyan)). 
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Figure 72. Salinity profiles for CTD station 65. Profiles shown are CTD initial (blue), 
and ocean model prediction seven days later (with AMPS forcing (black), 
observed forcing (red), AMPS wind speed increased by factor of two (green), and 
observed wind speed increased by factor of two (cyan)). 
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Figure 73. Potential density profiles for CTD station 65. Profiles shown are CTD initial 
(blue), and ocean model prediction seven days later (with AMPS forcing (black), 
observed forcing (red), AMPS wind speed increased by factor of two (green), and 
observed wind speed increased by factor of two (cyan)). 





3. CTD Station 92 


Figure 74 shows friction speed for the four different scenarios run at CTD station 
92. Observed forcing resulted in a higher friction speed than AMPS forcing. This trend 


is also shown when both AMPS and observed wind speeds were increased. 


A relatively large positive temperature flux area for the AMPS forecast scenario 
(Figure 75 a.) occurred at the end of day one into day two. For observed conditions 
(Figure 75 b.), large positive temperature flux areas occurred immediately on day one, 
between days three and four, and six and seven. For each of the increased forcing 
scenarios (Figure 75 c. and d.), the larger positive temperature flux areas occurred deeper 


in the water column, with very weak positive temperature flux near the surface. 
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For both AMPS scenarios (forecast and increased wind speed) the weaker forcing 
appeared to be unable to overcome the temperature flux. Figure 76 a. and c. show 
varying salinity flux (between positive and negative) throughout the water column. As 
before, buoyancy flux (Figure 77 a. and c.) followed salinity flux. The ocean model 
predicted ice thickness (Figures 78 a. and 79 a.) showed very little change for both of 


these scenarios. 


Stronger forcing associated with both observed scenarios led to negative salinity 
and buoyancy fluxes (Figures 76 b. and d., and 77 b. and d.) in the upper water column 
over the entire period. The ocean model predicted ice thickness for both of these cases 


(Figures 78 b. and 79 b.) show an almost linear increase in ice thickness. 


Potential temperature, salinity, and potential density profiles for are shown in 
Figures 80, 81, and 82 respectively. Initial CTD profiles and ocean model prediction at 
the conclusion of the run are as described for CTD station 1. As with CTD station 1, 


profiles agree well with contour and ice thickness plots. 


At CTD station 92, the observed forcing scenario resulted in the least stable water 
column at the end of the seven day model run. Potential density difference across the 
pycnocline (Figure 82) decreased from an initial value of 0.05 kg/m* to approximately 


0.01 kg/m*. Note that AMPS scenarios at this station caused very little change. 
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Figure 74. Friction speed (u«) contour plots from surface to 250 meters at CTD station 92 
for Julian days 233 to 240 a.) AMPS forecast, b.) Observed conditions, c.) AMPS 
forecast wind speed increased by factor of two, and d.) Observed wind speed 
increased by factor of two. 
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Figure 75. Temperature flux (<w'T'>) contour plots from surface to 250 meters at CTD 


station 92 for Julian days 233 to 240: a.) AMPS forecast, b.) Observed conditions, 
c.) AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two. 


99 


<w’'S'>, CTD 92, AMPS Input 10° <w’S'>, CTD 92, Observation Input 


0 
15 
- 1 
-100 
E 
150 
-200 , 






















Nn sy 
€ a 
Fs a 
2893 237 25933«234.=«8SSSC«8GSS*C«BT~SC«8BSSC«BD~=~«CO 
of 2005 Day of 2005 
<w’S'>, CTD 92, AMPS Wind Speed Increased by Factor of 2 <w'S'>, SID. 92, Soservenon Wind Speed Increased by Factor of 2 x10° 
0 0 
va — 
E E E 
-{ 2 f J 
~ e@ a 
“233 237 a 233 237 
sed of 2005 a of 2005 
Figure 76. Salinity flux (<w'S'>) contour plots from surface to 250 meters at CTD station 


92 for Julian days 233 to 240: a.) AMPS forecast, b.) Observed conditions, c.) 
AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two (Note white shaded contours are areas of 
relatively large negative values). 
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Figure 77. Buoyancy flux (<w'b'>) contour plots from surface to 250 meters at CTD 
station 92 for Julian days 233 to 240: a.) AMPS forecast, b.) Observed conditions, 
c.) AMPS forecast wind speed increased by factor of two, and d.) Observed wind 
speed increased by factor of two (Note white shaded contours are areas of 
relatively large negative values). 
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Figure 78. Ocean model predicted ice thickness at CTD station 92, Julian days 233 to 
240, for a.) AMPS forecast, and b.) Observed wind forcing conditions. 
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Figure 79. Ocean model predicted ice thickness at CTD station 92, Julian days 233 to 
240, for a.) AMPS forecast wind speed increased by factor of two, and b.) 
Observed wind speed increased by factor of two. 
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Figure 80. Potential temperature profiles for CTD station 92. Profiles shown are CTD 
initial (blue), and ocean model prediction seven days later (with AMPS forcing 
(black), observed forcing (red), AMPS wind speed increased by factor of two 
(green), and observed wind speed increased by factor of two (cyan)). 
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Figure 81. Salinity profiles for CTD station 92. Profiles shown are CTD initial (blue), 
and ocean model prediction seven days later (with AMPS forcing (black), 
observed forcing (red), AMPS wind speed increased by factor of two (green), and 
observed wind speed increased by factor of two (cyan)). 
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Figure 82. Potential density profiles for CTD station 92. Profiles shown are CTD initial 
(blue), and ocean model prediction seven days later (with AMPS forcing (black), 
observed forcing (red), AMPS wind speed increased by factor of two (green), and 
observed wind speed increased by factor of two (cyan)). 





E. OCEAN MODEL CONCLUSIONS 


Overall, it was found that observed forcing conditions compared to AMPS forcing 
conditions (whether actual observed or increased) led to greater change in the water 
column. The most dramatic example of this was at CTD station 92 (least stable water 
column), where weak AMPS forcing resulted in nearly no change in the stability, while 
observed conditions resulted in much more instability. It was also noted that at all three 
CTD stations, not one of the four scenarios was able to overcome the pycnocline and 


drive the water column to complete instability. 
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At CTD station 1, AMPS forecast and observed forcing was found to result in 
weak temperature fluxes in the water column, allowing for more ice growth. These 
values were 6 cm for AMPS forcing and 24 cm for observed forcing. The two increased 
forcing scenarios resulted in stronger temperature fluxes and decrease in ice thickness of 


6.5 cm for AMPS increased forcing, and 9.5 cm for observed increased forcing. 


At CTD station 65 wind forcing did not result in temperature fluxes strong 
enough to initiate ice melt. In all scenarios, the ocean model predicted ice thickness 
increases of 13 cm for AMPS forcing, 15 cm for AMPS increased forcing, 26 cm for 
observed forcing, and 28 cm for observed increased forcing. While the temperature 
fluxes did not cause ice melt, for increased forcing scenarios they did serve to limit ice 
growth, with increased forcing scenario ice thickness only 2 cm more than AMPS and 


observed. 


At CTD station 92, weak AMPS forcing led to weak temperature flux and little to 
no change in model predicted ice thickness. An increase of only 1 cm was seen for 
AMPS forcing and decrease of only 0.5 cm for AMPS increased forcing. Both observed 
scenarios resulted in stronger temperature fluxes, but not enough to promote ice melt. 
Model predicted ice thickness showed an increase of 15 cm for observed forcing and 13 


cm for observed increased forcing. 


Average change in potential density difference across the pycnocline at CTD 
stations 1, 65, and 92 were 0.075 kg/m° , 0.06 kg/m’, and 0.04 kg/m* respectively. As can 
be seen in Figure 25, wind speed decreased at each of these stations, with the weakest 
values at station 92. If forcing at CTD station 92 had been comparable to the other two 
CTD stations, there is a good possibility that compete instability would have been 


achieved resulting in overturning of the water column and deep convection. 
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IV. CONCLUSION, SUMMARY, AND RECOMMENDATIONS 


The purpose of this research was to explore the usefulness of NWP forecast data 
as an input into a simple ocean model, and the feasibility of using the ocean model output 


as a planning tool during reactive situations. 


The first step to accomplish this was to compare archived AMPS forecast model 
data with observed conditions during collected during the 2005 MaudNESS deployment 
to the Maud Rise area in the eastern Weddell Sea. Comparisons consisted of time series 


and statistical analysis. 


AMPS was found to have no skill in forecasting downwelling long wave radiation 
and relative humidity. It is believed that this was due to poor handling of cloud and 
moisture parameters. Setting relative humidity with respect to ice to saturation, and 
converting to relative humidity was found to produce slightly more accurate relative 
humidity values. It was also noted that AMPS consistently over-forecast near surface air 


temperature, which led to high values for specific humidity. 


AMPS surface pressure, near surface wind speed and direction, and downwelling 
short wave radiation all performed well. This led to the conclusion that, while there were 
problems with downwelling long wave radiation and relative humidity, AMPS seemed to 


well represent the overall synoptic situation. 


AMPS latent and sensible heat flux forecast skill is undetermined. This is due to 
the lack of observed surface temperature which led to the assumption (for bulk methods 
used to calculate observed heat fluxes) that surface temperature was equal to near surface 
air temperature. Comparison of AMPS and observed heat fluxes after this assumption 
resulted in no correlation between the two (with AMPS values appearing to be over- 
forecast). When AMPS bulk parameters (relative humidity, wind speed, air temperature, 
and surface pressure) were used to calculate heat fluxes using the bulk method, AMPS 


and observed values showed some correlation. The AMPS bulk method showed skill in 
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the latent heat forecast, despite no relative humidity skill, because of the dependence of 
latent heat flux on specific humidity difference and wind speed, which were predicted by 


AMPS with skill. 


The next step was to use AMPS forecasts and observed conditions as input into a 
simple one dimensional ocean model. This model was initialized with three different 
temperature and salinity profiles from CTD casts performed during MaudNESS. CTD 
station 1 was the farthest away from the Maud Rise and in the most stable water column. 
CTD station 65 was located over the Rise (in the Taylor column) in a marginally stable 
water column. The final station was CTD station 92, located to the southwest of the 


Maud Rise in the "warm pool’, and least stable water column. 


At every station, observed conditions (both actual and increased) resulted in the 
least stable ocean model output. The most dramatic example was at CTD station 92, 
where weak AMPS forcing resulted in almost no change in potential density profile, 
while observed conditions nearly resulted in eliminating the potential density difference 


across the pycnocline. 


One trend noticed, was the stabilizing effect of ice melt on the water column, 
particularly in the increased forcing scenarios at CTD stations 1 and 65 shown in Figures 
56 - 59 c. and d., 61, 65 - 68 c. and d., and 70. It was found that a strong wind event 
would begin mixing warm, relatively salty water into the mixed (or upper) layer, 
resulting in ice melt. The fresh water from the ice melt would then cause the upper layer 
to decrease in salinity and would curtail mixing, allowing freezing to resume and start the 


process Over. 


It was known that a simple one dimensional model could not accurately describe 
an area with such complicated variability as the Maud Rise region of the Weddell Sea, 
and was not the purpose of this study. As a sensitivity indicator, the ocean model proved 
extremely valuable. The utility of such an ocean model using NWP forecast model 


guidance for short term, reactive situations is entirely feasible, and should be explored for 
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mid latitude open oceans (for example to study effects on sound speed profiles and ASW 
applications). For long term situations, or where high accuracy is important, a more 


complex model would be required. 


It is recommended that further studies be conducted on modeled versus observed 
latent and sensible heat fluxes in polar regions with accurate observed surface 
temperatures. Also more research should be done on the utility of setting relative 
humidity with respect to ice to saturation in polar NWP models, which could result in a 
simpler, cheaper, and more accurate method of calculating relative humidity in a polar 


environment than current NWP model parameterizations provide 


In addition to using the Polar MM5, AMPS is also being run with the Weather 
Research and Forecasting Model (WRF). If the opportunity presents itself (e.g., 
additional expeditions to the Maud Rise/Weddell Sea area), similar research should be 


conducted comparing AMPS WRF to observed conditions. 
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